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Abstract 

Progress in manufacturing technology has allowed us to probe the behavior of 
devices on a smaller and faster scale than ever before. With increasing miniaturiza- 
tion, quantum effects come to dominate the transport properties of these devices, 
between collisions, carriers undergo ballistic motion under the influence of local 
electric and magnetic fields. The often surprising properties of quantum ballistic 
transport are currently elucidated in "clean" atomic physics experiments. From a 
theoretical viewpoint, the electron dynamics is governed by ballistic propagators 
and Green functions, intriguing quantities at the crossroads of classical and quan- 
tum mechanics. Here, we briefly describe the propagator method, some ballistic 
Green functions, and their application in a diverse range of problems in atomic 
and solid state physics, such as photodetachment, atom lasers, scanning tunneling 
microscopy, and the quantum Hall effect. 
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Figure 1: Timescales in semiconductors. 



1 Physics in small dimensions 

The laws of quantum mechanics provide a means for a successful interpretation of mea- 
surements and experiments on a microscopically small scale. It goes without saying 
that we can never "observe" directly what's going on in an atom or molecule. To un- 
derstand what nature is telling us we must learn its language. Its grammar follows the 
mathematical rules of quantum mechanics. Without mathematical tools we would not 
be able to describe intriguing processes such as, for example the mapping of quantum 
states of light to intrinsic atomic states. Fortunately, the necessary formalism is often 
powerful, elegant and quite easy to comprehend. In our contribution we will demon- 
strate the revival of an established mathematical tool in the important field of quantum 
technology. This tool is known under the name of Green function or Green's function 
in honor of George GreenQ 

Many problems in electrodynamics, hydrodynamics, heat conduction, acoustics, 
etc., require the solutions of inhomogeneous linear differential equations. It is there 
where Green functions come to full power. The corresponding mathematical approach 
is the same in all branches of physics — as long as we are dealing with linear, ordinary 
or partial differential equations. In quantum mechanics, Green functions enjoy the ad- 
vantage of having a physical meaning: The single-particle Green function G(r,x';E) 
is the relative probability amplitude for a particle to move with energy E from an ar- 
bitrary point r' to another point r. Probability amplitudes are known to be essential in 
all kind of quantum problems. In his book on 'The Character of Physical Law' Feyn- 
man |1 ] notes that ". . . everything that can be deduced from the ideas of the existence of 
quantum mechanical probability amplitudes, strange though they are, will work, . . . one 
hundred percent " 

In this tutorial we present some basic features of electron and atom motion in ex- 
ternal fields. External electric, magnetic and even gravitational fields are well suited to 
control the motion of particles in quantum devices. It is not our purpose to dive into 
the technical depths of ultrasmall electronics research and technology. We only want 

1 Green was an almost entirely self-taught English mathematician and physicist who in 1828 published an 
essay entitled: "On the Applications of Mathematical Analysis to the Theories of Electricity and Magnetism." 
In this essay he obtained integral representations for the solutions of problems connected with the Laplace 
operator. 
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to illustrate how useful single-particle Green functions and propagators can be for ba- 
sic problems in quantum technology whenever a microscopic description of quantum 
transport is necessary. A single particle description is appropriate for devices with low 
particle densities where interaction processes can be neglected. But they are also useful 
in more general cases. Indeed, microscopically small particles travel freely on length 
scales of the order of the free mean path I, which is the distance that an electron travels 
before its initial phase is destroyed for whatever reasons J2JE1. The mean free paths 
depend strongly on the material under consideration and they are much affected by 
temperature. Particles that travel freely are called ballistic particles. 

In Fig.[TJwe show timescales for the motion of electrons in a typical semiconductor. 
We invoke the uncertainty principle Ap Ax ~ h to get a feeling for the time domain of 
ballistic motion. Motion with well defined momentum requires Ap <C p- For a quasi- 
classical description, the electrons are required to be well localized compared to the 
mean free path, i. e., Ax -C £ = px/m with x being the time for ballistic motion (see 
Fig.l). It follows that ApAx <^ipi — p 2 x/m. Suppose we have thermal electrons; then 
we can use the equipartition law of classical statistics, p 2 x/m ~ kgTx with kg being the 
Boltzmann constant and T the temperature. Comparing with the uncertainty principle 
we obtain x ^> h/{ksT) which for room temperature is of the order 10~ 13 to 10~ 14 
seconds. For thermal electrons this last inequality is frequently met. However, for 
electrons moving with high energies of the order e V the inequality is no longer fulfilled. 
In this case one has to treat ballistic transport fully quantum mechanically. In the 
following we present examples where quantum transport is essential. First, however, 
we review some useful mathematical tools in a nut shell. 

Clearly, we cannot cite all relevant literature in this field. This would be an impos- 
sibly difficult and lengthy job. But the interested reader will find a wealth of literature 
in the research articles cited in this tutorial. 

2 Propagators and Green functions 

A time-dependent treatment of the flow of charge carriers is based on the time-dependent 
Schrodinger equation. In this context it is useful to summarize a few aspects of the 
initial-value problem for a wave function known at t = to, 

\lf(r,t = t ) = Vo(r) = {r\y{t = t )}. (1) 

The corresponding ket vector evolves according to the basic law of quantum 

mechanics, 

(ihd,-H)\ W (t))=0. (2) 

The formal solution of (O is conveniently written in terms of the time evolution oper- 
ator U(t,to): 

\w(t)) = U(t,t )\Y(to)). (3) 

In coordinate space (O reads 

<r|y(/)) = / dV(r|t/( f , fo )|r')(rV(fo)). (4) 
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The integral kernel of (HJi is called propagator K: 

K(r,t\r',t ) = (r\U(t,t )\r'). (5) 

Obviously, K(r,t |r',fo) is the time evolution matrix f/(f ,?q) in coordinate space repre- 
sentation. From the last two equations, we have 

limtf(r,f|r',f H<5 3 (r-r'). (6) 

l->'0 

The (time-) retarded Green function must vanish for t < fo. It is usually defined by 

G(r,?;r',/ ) = -e(t-t )K{r,t\r',t ) . (7) 
m 

From this definition, and the fact that the propagator is a solution of the time-dependent 
Schrodinger equation, the retarded Green function is seen to satisfy the differential 
equation 

(ihd t -H)G(r,t;r',t ) = <5 3 (r-r')<5(f-f ) (8) 

The delta function 8(t — fo) on the right hand side of © originates from the step func- 
tion ©(f — fo) in the definition of G. 

For time-dependent Hamiltonians, the propagator will depend separately on f and 
fo. For time-independent Hamiltonians, the propagator depends only on the time dif- 
ference f — fo. In the latter case the Laplace transform of the propagator 

G(r,r';£) = - I dtexp(iEt/h)K(r,t\r',0) (9) 
in Jo 

generates the energy (-dependent) Green function, 

1 



G(r,r';£) = lkn 

7)^0+ 



E-H + irj 



(10) 



G(r,r';£) is the amplitude for travel of a particle from r to r' out of a point source 
and, as a function of energy. This feature will emerge if we evaluate ( TTOb explicitly. We 
should also mention that the appearance of the infinitesimally small, positive imaginary 
term irj in ( TTOb has a simple reason: To enforce convergence of the integral in (O one 
has to replace E by E + irj . The physical meaning of such a small shift into the complex 
energy plane becomes evident if one evaluates (TToT > for a free particle. The result (4), 

G fr ee(r,r,£)=- i -^ j—yj , (11) 

is well known from scattering theory: For r' fixed and r variable, Gf ree describes an 
outgoing spherical wave that originates from a point source at r = r'. Had we Fourier 
transformed the time-advanced Green function instead of the time-retarded Green func- 
tion, we would, of course, have ended up with an incoming spherical wave instead of 
an outgoing wave. 

Propagators contain all necessary information about the motion of a particle. Un- 
fortunately it is not always possible to find a closed-form solution for K or G. For 
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Figure 2: The Moshinsky shutter: The shutter is removed at fo = 0. The quantum 
particles start propagating towards the screen. When will the particles arrive at the 
screen? The distribution of the traveling particles is shown on the right hand side where 
the full line corresponds to the quantum result (17) and, where the classical propagation 
is represented by dotted lines. 



potentials which are at most quadratic in the coordinates, the propagator assumes the 
canonical form E) 



*(r,r|r',0) =A(t) exp [iS c i(r,r';t)/h] , 



(12) 



where S c i is the corresponding classical action, and where A(t) is a time-dependent 
factor independent of the particle's position. However, nonquadratic potentials such as 
the Coulomb potential generally do not have the canonical form (fT2l . Explicit expres- 
sions for propagators can be found, for example, in ll5l 161171 l8l l9l [TOl fTTI . 



The Moshinsky shutter: This example illustrates how the free propagator 



tfftee(r,f|r',/o) = 



1 1/2 



2nih(t-t ) 



exp 



im(r — r ) 



'\2l 



2h(t-t Q ) 



(13) 



is used to solve an initial value problem that describes the flow of quantum particles. It 
is of interest in the context of the quantum mechanical propagation of a signal. Moshin- 
sky lfl2l has analyzed the spreading of such a signal. He considered a monochromatic 
beam of noninteracting particles of mass m and energy E^. The particles are supposed 
to move parallel to the x-axis from left to right. The beam is stopped (and absorbed) 
by a shutter at x — (see Fig. |2j. The signal is given at t = to = when the shutter is 
opened. 

The sudden removal of the shutter marks the beginning of a "quantum race" where 
the particles run along the positive x-axis. In order to elucidate the spreading of the 
signal all one has to do is to calculate \jf(x,t) starting with 



y/ (x) = yr{x,t = 0) = 0(-x)e 



ik.x 



(14) 
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Using one finds 



(x|y/(f>0)) =M{x;k;ht/m). 



(15) 



where the Moshinsky function M is defined in terms of the complementary error func- 
tion 03), 



M(x;k;z) = ^exp{ikx-ik 2 z/2) erfc 



x — kx 



(2iT)V2 



(16) 



with i 1 / 2 = exp(i^/4). 

An interesting property of <JT5j is revealed when we evaluate the particle number 
probability. Introducing u = {Tikt / m — x) / {nhl / m) 1 1 2 , we obtain 



(x\ W (t))\ 



+ C{u) 



■S(«) 



(17) 



The functions C(«) and S(m) are the well-known Fresnel integrals 01311 . The corre- 
sponding probability pattern is called diffraction in time because it arises when the 
shutter is opened for a finite time t . Although transient effects are important by them- 
selves 1T4) we won't discuss them in more detail here. In what follows, we will discuss 
stationary quantum transport. 



3 Quantum sources 

In real-space representation, propagators and Green functions describe the motion of 
quantum particles from some initial point r' to a final point r. But where do the parti- 
cles come from? One may think of two different situations: i) the particles have been 
around all the time like electrons in an atom, or, ii) the particles are generated by a 
source, a situation which is quite familiar from scattering theory where a beam of par- 
ticles is generated by an accelerator in a region far away from the target. In mesoscopic 
physics and nanotechnology, however, there is usually no such large spatial separation. 
Let us motivate the introduction of coherent quantum sources of particles and illustrate 
their properties by means of an example. 

3.1 Photoelectrons emitted from a quantum source 

We may consider the photoeffect as a two-step process as illustrated in Fig. [3] The 
time evolution of the emitted electron is of course governed by the rules of quantum 
mechanics. In the absence of any interaction between photon and electron, the electron 
under consideration is attached to the atom and is described by the bound-state wave 
function y/atom(r) . Let us consider a dilute gas of independent atoms where the interac- 
tion of the photoelectron with neighboring atoms can be neglected. In the presence of 
a photon field this wave function will obtain a small scattering component y/ S c(r) that 
allows the electron to leave the atom. For a dipole-allowed transition, the dipole oper- 
ator D(y) oc £ ■ r is responsible for transferring the electron from its initial bound state 
IVatom) to a continuum state |i// sc ) . Under steady-state conditions with many atoms 
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Figure 3: Two steps to create a photoelectron: Step 1 (left panel): The photon trans- 
fers its energy to the initially bound electron. Step 2 (right panel): The photoelectron 
escapes from the absorption region. 



(each having the Hamiltonian //atom) and weak monochromatic light we must solve the 
problem: 



E - //atom - Z/rad - -D • (fl + fl 1 " 



ID 



lv sc )|o) =° 



(18) 



where E = £atom + hv is the energy sum of electron and photon. The unperturbed 
Hamilton operator of the radiation field with field operators a and cfi is denoted by // ra d, 
with the zero-point energy being subtracted. As usual, 1 1) characterizes the presence of 
the photon and |0) its absence after absorption. Projection onto the zero-photon state 
(0| yields the desired equation for the scattering solution, 



£-#atom(r) y/sc(r) = /)(r)v/atom(r) ee a{r). 



(19) 



We can interpret D\y/) as a source function \a) for the photoelectrons: The dipole op- 
erator prepares the electron in a continuum state but can be neglected once the electron 
has left the atom. 



3.2 Currents generated by quantum sources 

The last two equations can be generalized to a situation where the scattered particle 
experiences some final-state interaction. For example, the presence of a final-state 
Coulomb interaction or of an external field can be readily taken into account in Eq. ( fT9l 
by writing 

£-// -W(r)]v/ sc (r) = c7(r), (20) 

where Ho is the Hamiltonian of a free particle and where W(r) represents the interac- 
tion of the emitted particle (for example the photoelectron) with its environment. Here 
and in the following we omit the hat symbol for the operator (Ho = Ho). In analogy to 
other inhomogeneous field equations, e. g., Maxwell's equations, the right-hand term 
<r(r) in d20t is again identified as a source for the scattered wave y/ sc (r). 

We now turn to the mathematical aspects of d20l ). Introducing the energy Green 
function G(r,r / ;£) for the Hamiltonian H defined via fl31 : 

[E -H - W(r)} G(r,r';£) = S 3 (r -r') , (21) 
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a solution to (120b in terms of a convolution integral reads: 



Yfc(r) = |dVG(r,r';£)c7(r') 



(22) 



In general, this result is not unique. However, any two solutions Ysc( r ) an d Ysc( r ) 
differ only by an eigenfunction y/hom( r ) of the homogeneous Schrodinger equation, 
with H = Hq + W and H\j/hom{ r ) = E\f/^ om (r). The ambiguity in y/ sc (r) is resolved 
by the demand that G(r,r';E) presents a retarded solution characterized by outgoing- 
wave behavior as r — > °°. Formally, this enforces the same choice as in Eq. ( TlOb . It is 
then easy to decompose the Green function (TlOb into real and imaginary parts 



G(r,r';£) 



PP 



1 



E-H 



m8(E-H) 



(23) 



where PP(. . .) denotes the Cauchy principal value of the energy integration. 

Defining the current density in the scattered wave in the usual fashion by j(r) = 
^3[Vsc(r)* Vy/ sc (r)]/M (where for simplicity we omitted the vector potential A(r), see 
lfl6l ). the inhomogeneous Schrodinger equation d20T > gives rise to a modified equation 
of continuity lfT71[T8l : 

V.j(r) = -^3[cr(r)> sc (r)] , (24) 

where 3 [x] stands for the imaginary part of x. Thus, the inhomogeneity (j(r) acts 
as a source for the particle current j(r). By integration over the source volume, and 
inserting d22b . we obtain a bilinear expression for the total particle current J(E), i. e., 
the total scattering rate: 



d 3 r / d 3 r'(j(r)*G(r,r';£)(j(r') 



(25) 



Some important identities concerning the total current J(E) are most easily recognized 
in a formal Dirac bra-ket representation. In view of d23l , we may express J(E) by 



2 In 
J{E) = --3 [(a \G\ a)} = —{a \8{E-H)\a) , 

from which the sum rule immediately follows IT8l : 

jyEJ(E) = *£(o\o) = ^J d 3 r|a(r)| 2 , 
(provided this integral exists). 



(26) 



(27) 



3.3 Recovering Fermi's golden rule 

In order to connect Eq. d25l l to the findings of conventional scattering theory, we display 
J(E) in an entirely different, yet wholly equivalent fashion. Employing a complete 
orthonormal set of eigenfunctions | i//fi) of the Hamiltonian H, 8 (E —H) \ i//fj) = 8 (E — 



2 mm 



Photon energy hv 



Figure 4: Near-threshold detachment of oxygen ions: O -^0 + e~ in the presence of 
a homogeneous electric force field F = eS . The two possible classical trajectories for a 
photoelectron leading from the source (marked by S) to any destination will give rise to 
interference on a distant detector screen. The fringe pattern in the current distribution 
depends sensitively on the energy. By counting the number of fringes the binding 
energy Eq of the outer electron can be determined from Einstein's law l22l l23l l24l . 

£fi)|l//fi) follows, and replacing \a) = ZD(r) | V^atom) dHJl, we may formally decompose 
d26l i into a sum over eigenfunctions: 

J(E) = —Y, S ( E - E ^)\(¥n\D{r)\xi/ atom )\ 2 ■ (28) 
" fi 

Thus, Fermi's golden rule is recovered. Another noteworthy consequence of dZSb and 
d26l i emerges in the limit of pointlike sources, a(r) ~ C8(r — R). We then find ifTTl 

J(E) = --|C| 2 3[G(R,R;£)] = ^|C| 2 «(R;E) , (29) 

where n(R;E) = £ fi 8(E — £fl)| yAfl(R) | 2 is the local density of states of H at the source 
position R. Equation ( 1291 forms the theoretical basis of the Tersoff-Hamann descrip- 
tion of scanning tunneling microscopy |fT71[T9l . The advantage of the formulation in 
terms of quantum sources over the traditional Fermi's golden rule approach (which in- 
volves an integral over the final states) is that it emphasizes the dynamical aspects of 
the propagation in real space and opens the possibility to a semiclassical calculation of 
photocurrents with closed-orbit theories 1201I2TI . 

3.4 Photodetachment and Wigner's threshold laws 

To find out how we can use the formalism for real physics we continue our discussion of 
the photoeffect. Applying the photoelectric effect to negative ions means that the emit- 
ted electron only weakly interacts with the remaining neutral atom. Just as in Young's 
double-slit experiment, the fringe pattern in the current profile can be interpreted as 
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interference between the two classical trajectories, here of a particle in a constant force 
field l24l l25l . From the interference pattern one can determine the kinetic energy of 
the electrons and plot it against the photon energy to check Einstein's law (right panel 
of Fig. Near threshold, the photoelectron has very little kinetic energy and, hence, 
a large de Broglie wavelength. In the absence of a final-state Coulomb interaction the 
relevant Green function is that of a particle falling freely in a constant field ll24l l26l . 
Experimental results are reported in Ref. 1231 . Fig. 4. They show a highly accurate 
verification of Einstein's law which can be used to obtain the binding energy of O 
with unprecedented accuracy. 

As noted above, we can interpret D(r) | y a tom) as a source | cr) for the photoelectrons. 
Expanding the source in terms of multipoles <7/„,(r) and, taking into account that for 
O the photoelectron leaves the atom near threshold in an s-wave continuum state, 
we retain only the / = component of the source by writing c(r) = Coot/) / V47T . In 
the absence of external fields, the free Green function (TTTb . a outgoing spherical wave, 
yields after multipole expansion 



J(E) °c k 



rdr 

— sm(kr)GQo(r) 
o k 



(30) 



For E = h 2 k 2 /(2m) — ^ it follows that J(E) k. This result is independent of the 
form of the atomic source and it reflects Wigner's threshold law Il24ll27l . It provides 
an alternative way to determine the electron affinity of a negative ion. 

A similar analysis applies to the more complicated Green function in an electric 
force field F = e<§, one of the few quantum problems in more than one dimension that 
have exact solutions: 

G(r,o;£) = [Ci(a+)Ai'(a_) -Ci'(a+) Ai(a_)l . (31) 

Here, j3 3 = m/(2/iF) 2 , a± = -/3[2£+F(z±r)], and Ai(u) and Ci(u) = Bi(it) +i Ai(«) 
denote Airy functions fl3l . It leads to a modified Wigner law for the s-wave absorption 
cross section near threshold lfl8ll : 

J{E) °c Ai'(-2j3£) 2 + 2/3 £ Ai(-2j3£) 2 . (32) 

A static electric field opens up a sub-threshold (E < 0) tunneling regime that has been 
confirmed by experiment l28l . 

Emission of particles from pointlike sources has been considered in the literature 
||29l long before the advent of mesoscopic physics. It was Schwinger [ 30 1 who intro- 
duced sources as a means of describing quantum dynamics in the context of emission 
and absorption of light. 



4 Spatially extended sources: The atom laser 

Atomic electron sources are usually sufficiently small to be considered pointlike. A dif- 
ferent situation arises when particles are coherently emitted from an extended region 
in space. An example for such a "fuzzy" source is the continuous atom laser, a beam 
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Figure 5: Left panel: A Gaussian source of freely falling particles can be replaced by a 
virtual point source of particles with the same energy, located upstream from the actual 
extended source. Right panels: Size dependence of the beam profile for Gaussian BEC 
sources with different widths a. 

of ultracold atoms fed by a Bose-Einstein condensate (BEC) 1311 . In the experiment, 
only atoms in a specific Zeeman substate (to = —1) are magnetically trapped and form 
a BEC. Application of a suitably tuned radiofrequency (RF) field will cause transitions 
into another magnetic substate of the atoms (m = 0) that is not subject to the trapping 
potential. Under the influence of gravity, these "outcoupled" atoms fall freely from 
the trap region and form a coherent, continuous atom laser "beam." In our language, 
the macroscopic BEC wave function y/b(r) serves as the source cr(r) and corresponds 
to the atomic bound state V / atom(r) in Eq. (fT8l l. whereas the outcoupled beam y/(r) of 
accelerating atoms takes over the role of the scattered wave y/ sc (r), akin to photode- 
tachment in an electric field (Sec. l3.41 >. From a theoretical viewpoint, the only essential 
difference is the macroscopic size of the source. 

For ideal, non-interacting atoms, the ultracold cloud populates the ground state of 
the nearly parabolic trapping potential, leading to a Gaussian density profile in the 
BEC. For simplicity, we assume an isotropic distribution: 

ff(r) = Mly (r) = MW exp(-r /(2a 2 )). (33) 

Here, hQ. denotes the strength of the transition-inducing oscillating RF field. The pa- 
rameter a describes the width of the source (which is related to the field gradient in the 
trap), and No = a~ 3,/2 7r~ 3 / 4 denotes the proper normalization from the condition 

|d 3 r|v/o(r)| 2 = l. (34) 

To obtain expressions for the currents generated by a Gaussian source, we work in the 
time-dependent propagator representation (see Eq. (0). The beam wave function i//(r) 
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Figure 6: Left panels: Atom-laser beam profile from a Rb BEC of size a = 0.8 |im 
at different detuning energies Av = E/h A3 61 . The first series shows the beam profile 
for non-interacting particles, whereas the next series includes interactions due to 500 
atoms, which lead to a transverse substructure l33l . Right panel: Anisotropic trapping 
frequencies cause a strong modulation of the total particle current J(E). 



then may be written 

y(r) = -iiW jT 'dte iE, l h j dV tffieid Mr' ,0)e- /2/(2a2) . (35) 

It is possible to carry out the integration over the source volume. With negligible 
corrections, outside the source the integral d35l l assumes the form 

W (r) =Si2(2^ fl ) 3/2 e-™ 2£ / fi2+m2F2 " 6 /( 3fi4 )G(r,-g fl 4 ;£), (36) 

where G(r,r';£') denotes the energy Green function for uniformly accelerated parti- 
cles (f3TT >- This expression displays a remarkable feature of the beam wave function 
y/(r) originating from a Gaussian source: The extended source can be formally re- 
placed by a virtual point source of the same energy, albeit at a location shifted by 
r' = — mFa 4 / (2h 2 ) from the center of the Gaussian distribution (see Fig. [5]). Expres- 
sions for the beam profile and currents are then conveniently found from the analogous 
expressions for a point source by performing the indicated shifts. 

As an immediate, and somewhat surprising, consequence of the concept of a virtual 
source, the beam profile shows a sharp fringe pattern that results from the interference 
between the two virtual paths in Fig. [5] The number of fringes depends sensitively on 
the size a of the source, as displayed in Fig. [5] In the limit of extended Gaussian sources 
with E < mF 2 a 4 / (2h 2 ), the virtual source turns into a tunneling source (as discussed 
in greater detail in the following section), and the beam profile itself becomes Gaussian 
11261 . The spectrum of the total particle current J(E) as a function of the detuning of 
the RF field E = hAv then may be written in the suggestive form 

J(E)^^Jd\\a(r)\ 2 8(E + Fz), (37) 
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Figure 7: Beam profile for simultaneous output coupling with two different radio 
frequencies. The outcoupling frequencies Avi,2 —E/h are (a) ±0.5 kHz, (b) ±1.0 kHz, 
and (c) ±2.0 kHz. The number of longitudinal interference fringes is proportional to 
the difference in the detuning frequencies. Parameter: a = 0.8 jim, F = m^og, with 
g = 9.81 m/s 2 , and WRb = 87 u. 



an expression that has a simple geometrical interpretation. For extended sources, the 
energy dependence of J(E) reflects the source structure: By the resonance condition 
E + Fz = 0, the total current probes the density | \f/Q (r) | of the BEC on different slices 
across the source. Finally, we note that the approximation ( T37T ) obeys the sum rule d27T i 
for the total current J(E). 

In an actual atomic BEC, the repulsive interactions between atoms lead to a broad- 
ening of the condensate. For most cases, the inclusion of the interactions via a mean- 
field approach is sufficient. The repulsive forces of the much denser BEC act on the 
outcoupled atom beam and lead to a further splitting of the beam profile, as shown in 
Fig- 151 Also the total current is modified by the interactions J33). Both effects have 
been observed experimentally. Non-isotropic trapping frequencies and currents from 
higher trapping modes allow to control the shape and rate of the atom laser 1331 . 

Fig.|7]shows the "dripping quantum faucet", which is produced by superposition of 
two laser beams with slightly different energy that are outcoupled from the same BEC 
[32]. It is not surprising to see that rotating BECs which sustain vortices are described 
in terms of rotating Gaussian sources with nodal structures 11241 . 



The quantum theory of scattering is not limited to asymptotic problems where particles 
are generated (and observed) far away from the scattering region. A prominent can- 
didate for scattering at finite distances is the Scanning Tunneling Microscope (STM). 



5 Ballistic tunneling: STM 
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Figure 8: Electric field emission out of (or into) the sharp apex of a conducting wire. 
The bright spot symbolizes the Gaussian density profile along the central tunneling 
path between tip and surface of the sample. 

There, an electric current flows down a macroscopic wire that ends in a sharp tip. Its 
apex can be viewed as a source of electrons which leave the tip by tunneling due to 
the applied electric field between tip and sample surface. In some cases, the apex of 
the tip is ultra-sharp, consisting of a single atom. In an experiment the tip is slowly 
moved across the surface. In the constant current mode the tip is raised and lowered so 
as to keep the current constant. The raising and lowering process produces a computer- 
generated contour map of the surface 1371 [38l . The method is capable of resolving 
individual atoms and works best with conducting materials. Electrons drawn from the 
apex of the tip (see Pig- [SJ exhibit dynamically forbidden motion because the electron 
transfer between tip and surface occurs via field-driven tunneling, confining the current 
to a narrow filament with Gaussian profile that samples the surface. 

It is straightforward to model the apex of an STM tip as a source (or sink) of elec- 
trons. To be specific, we consider here a conducting sample surface that harbors a 
two-dimensional electron gas; in practice, the band of surface states on the densely 
packed, smooth Cu(lll) surface has been exploited for this purpose 11351 1361 . The 
STM tip will emit a spreading surface electron wave that is scattered at adsorbed sur- 
face atoms (adatoms). Since the electrons are slow, s-wave scattering prevails that 
can be modelled by a short-range potential. In this case, one has an analytic solution 
for the scattering problem which forms the basis for the calculation of the corrugation 
(surface roughness). As a result of the analytic approach we will show that scattering 
resonances play an essential role for resolving atoms and detecting electron surface 
states. 

The current flowing through the STM tip is proportional to the local density of 
surface states, and therefore the imaginary part of the Green function at the tip posi- 
tion, 3[G(R,R;E)] ( 1291 . A route that leads conveniently to the Green function in this 
problem consists of the following three steps: 
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Figure 9: Sketch of the bouncing ball problem. An electron exits the three-dimensional 
tunnel at z = with energy E = 0. It then undergoes multiple reflections between 
the exit of the tunnel and the surface barrier at z = Zw before travelling into the solid 
(z > z w )- Adatoms (not drawn in the figure) are necessary to observe the resonance- 
induced ripples (see Fig. ITTb with energy E±. 

Step 1. One-dimensional problem: We first calculate the one-dimensional Green 
function G"(z,z';E) that corresponds to the model potential of Fig. [9] In this one- 
dimensional problem, the electron is allowed to tunnel in direction of the electric field 
(z-direction). The energy Eq corresponds to the bound (and unoccupied) surface state 
of Cu(l 1 1). The energy E of the tunneling electron is taken to lie in a band gap of the 
substrate (solid). As a result the electron faces a potential barrier at z = Z w an d bounces 
back and forth between tunnel exit (z = 0) and barrier. However, the electron can move 
freely with energy E± = E — Eq in the surface plane orthogonal to z = Z w - Because 
of inelastic scattering with phonons the electron will finally disappear in the solid. By 
assuming a point source at z = z! the Green function that belongs to the problem of 
Fig.|9]can be solved analytically in terms of Airy functions ll34l . 

Step 2. Three-dimensional background Green function: Since the STM is a three- 
dimensional device we must calculate the Green function in three dimensions. The 
uncertainty principle for momentum and position applied in the lateral (x—y) direction 
results in a tunneling spot (see Fig. [8) of finite width, and approximately Gaussian 
profile. The three-dimensional Green function G sym (r,r',z) for a particle moving in 
the potential U(r) = U(z) of Fig. [9] is obtained from its one-dimensional counterpart 
by integrating G" over all momenta hk± vertical to the field direction (i. e., parallel to 
the surface), 

G sym (r , r' , E) = G sym (z,z' , Ap,E) = ^ jT ' dkk ± J (k ± Ap)Gf (z,z' ,E - ^ ) , 

(38) 

with ./o(' ' ) being the usual cylindrical Bessel function of degree zero. The lateral 
distance between r and r' is given by Ap 2 = \p — p'\ 2 = (x — x') 2 + (y — y 1 ) 2 - Obviously 
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Figure 10: Schematic plot of tip (at z = z')> adsorbed atoms (at z = 0), and a strongly 
reflecting, clean surface (at z = z w ). 




Figure 11: Interference of waves scattered by corral atoms: Model calculation [34| of 
a corrugation plot at constant electric current. The structure of the circular ripples both 
inside and outside the corral can be related to the quantum bounce problem illustrated 
in Fig. [9] and discussed in the text. For comparison with experimental results obtained 
by Eigler's group [35 36] we show a similar setup with a quantum corral consisting 
of 48 iron atoms on a circle with radius 71.3 A adsorbed on a Cu(lll) surface (left 
panel). The corrugation of the adatoms is approximately 0.5 A and corresponds to a 
conductivity of Ob = 2.7 • 10~ 8 A/V. Note that there is no adsorbed atom in the center 
of the corral. 
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G sym (r',r / ,£') is independent of the lateral position (V,/) of the tip. Hence, for z' = 
const, G S ym represents a constant background corrugation. The full solution for the 
Green function with the adatoms present, is then obtained from 

Step 3. Dyson equation for G(r,r',£): We must now take into account the adsorbed 
atoms (see Fig. ITOb . Using the appropriate Dyson equation, we obtain an algebraic 
equation for the full Green function, 

n 

G(r,r',£) = G sym (r,r' ! £)+ £ G sym (r,r j ,E)(T(E)) jk G sym (r k y i E), (39) 

where the sum runs over all adatoms, and the r -matrix describing the effects of the 
multiple scattering events between the adsorbed atoms can be expressed using the 
background Green function G sym (r ; -,r^,£'). Details of the calculation of the Green 
function and the experimentally observable tunneling current J(R;E) d29l can be found 
in Ref. 11341 . A zero-temperature plot obtained from such a calculation (which typically 
takes a few minutes on a personal computer) is shown in Fig.Qj] We should point out 
that since G sym gives rise only to a uniform background current, the observed roughness 
of the surface is entirely contained in the r -matrix. 

6 Electrons in Electric and Magnetic Fields: 
The Quantum Hall Effect 

In this section we explore the strange and fascinating ways of electrons in electric and 
magnetic fields. Of particular importance here is the Hall configuration, where the 
electrons are confined to an effectively two-dimensional conductor in the presence of 
orthogonal electric and magnetic fields. The Hall geometry is displayed in Fig. [T2] 
Fig. Q~3] shows some of the classical paths followed by the electrons in the conduct- 
ing plane. Notwithstanding the complicated pattern of motion, all trajectories share 
the same distinctive behavior, uniform drift motion perpendicular to both fields with a 
characteristic velocity \'o = § 1 3$. 

Because of the universal drift motion, the current density in the Hall bar is simply 
proportional to the local density of states (LDOS) n(E) in the material, which in turn 
is related to the Green function in the corresponding external potentials (see Sec. 13.4b . 
Hence, from a mathematical point of view we are interested in finding the energy- 
dependent density of states for the moving electrons. For a purely magnetic field, 
the two-dimensional LDOS has a spike-like structure PT1 l42l . formally written as a 
superposition of discrete 5-distributions positioned at the Landau levels at E = (2k + 
1 )h(£>L, where col — eSSj (2m) denotes the Larmor frequency. 

n% D) (E) = ^-Y.8(E-ti03 L [2k+l]). (40) 

The addition of an electric field leads to important changes in the density of states 
for a purely magnetic field. As we know from Eq. d29l . the local density of states n(E) 



17 



Magnetic Field B 




Figure 12: Schematic view of a Hall bar. A current J x is flowing through a two- 
dimensional electron gas (2DEG) in the .Y-y-plane, which is orientated perpendicular 
to an external magnetic field SB. The deflected electrons at the sample edges produce a 
Hall voltage U y over the sample width W, which is measured along with the longitudi- 
nal voltage drop U x . 



is always linked to the imaginary part of the energy-dependent retarded Green function 
G(r = o,r' =o;E), 

n(E) = --3[G(o,o;E)\, (41) 

71 

which in turn can be expressed as the Laplace transform of the quantum propagator 
K{o,t\o,G) ©: 

G(o,o;£) = -/ &tz m, l n K{o,t\o,Q). (42) 
in Jo 

For r = 0, the two-dimensional quantum propagator for the Hamiltonian of a (spinless) 
electron in crossed fields, 



H 



(2D) 



Px + P? 1 



-m 



col (x 2 +y 2 ) - r_L ■ F_l - Pyxa> L + p x yco L , 



s * m 2m • 2 
where Fx = — e£± denotes the electric force in the jt-y-plane reads IfTTII . 

(2D) , \ imOOr [ iFh , , . 

K^zUo,t\o,0) = - - - . . — -exp<^ — ^ [co L t cot (co L t) - 1 



(43) 



(44) 



To account for the effects of the electron spin, we note that its interaction with the mag- 
netic field merely adds a spatially constant term to the Hamiltonian, and thus shifts the 
effective energy of the two spin populations by a fixed amount ijg/igj? = ±jgh(OL 
BP . The spin-dependent densities of states become 



(45) 



m \(E) =n I E± -gh<x>L 
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Figure 13: Drift motion in crossed fields: Trajectories drift with velocity \'o = £ j 33 
perpendicular to the fields. The motion of the electrons in two dimensions is a su- 
perposition of cyclotron motion with drift motion, resulting in trochoidal (cycloidal) 
trajectories. 



and the total LDOS including spin can be mapped back to the LDOS without spin: 
n^i(E) = n^(E) +n^{E). Thus, it suffices to evaluate n(E). For the discussion of 
the density of states in the Hall configuration it is useful to replace the trigonometric 
functions in the propagator (l44l by a sum. This can be done using the identity 

e X p[-acoth( z )] = 2e _ a - L (0) (2a)e - 2 ^ +1/2); 6) 
sinh(z) 4 fo 

which follows directly from the generating function of the Laguerre polynomials l!? (z) 
lfl3l . The Laplace transform d42T > then can be performed analytically, and we finally 
find for the density of states 



"<*" iE) = 2^t^°~^ mEk ™ 2 ' (4?) 
where the level width parameter 



T = F ± l (48) 



is related to the magnetic length / = ^Jh/{e38). H^(z) denotes the kxh Hermite poly- 
nomial [13] and is the effective energy shift for the kth Landau level 

E k =E-T 2 /{4h(D L )-{2k+\)h(0 L . (49) 

Interestingly, the density of states can again be interpreted as a sum over Landau levels. 
However, they now appear spread in energy, with a distribution that is isomorphic to 
the probability density of the corresponding eigenstate of a one-dimensional harmonic 
oscillator 

k(£)| 2 = ^^e-« 2 [H fc (§)] 2 (50) 

The total contribution of the A:th Landau level integrated over energy-space is readily 
available from the normalization of the oscillator eigenstates: 

'l/'. 'U ' /:/•• < / " ^ l M ^)| 2 = |§- (51) 
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Figure 14: The quantum Hall effect. Left panel: Schematic sketch of the experimental 
results adapted from |39] for p xy in GaAs/AlGaAs heterojunctions at T = 50 mK. The 
theoretical simulations 11401 (right panel) are based on the density of states d47l >. The 
dotted straight line is the classical prediction. 



This result reflects the quantization of each Landau level in a purely magnetic field 

At low temperatures, only the occupied states with energies smaller than the Fermi 
energy Ep of the system will contribute to the Hall current. Thus, we expect that 
the current density will be proportional to the integrated density of states N(Ep) with 
energies E <Ep, 

N(Ep) = J F 6En(E). (52) 

Within the Fermi gas model for a dilute gas of electrons in two dimensions it is straight- 
forward to obtain resistivity plots which bear remarkable resemblance to actual quan- 
tum Hall data l40l . An example is shown in Fig. [14] 

We finally mention the implications of our Green function model of the integer 
quantum Hall effect for the observed breakdown of the quantized conductivity at higher 
electric fields 11431 . Kawaji l44l 1431 14"6l studied the width of the quantized resistivity 
"plateaus" as a function of the electric current and thus the Hall field in the sample. 
He finds a characteristic power law for the shrinking of the plateaus, which can be 
expressed in terms of a critical electric field: 

4rit-Z? 3/2 . (53) 

The non-perturbative inclusion of the electric field via the Green function formalism as 
presented here leads automatically to this power law (see l40l[47l[48ll49l ) 

7 The semiclassical method 

It is usually rewarding to analyze quantum problems with the tools of classical physics. 
First, the semiclassical solution facilitates in many cases the understanding of the quan- 
tum solution properties. Second, semiclassics leads to approximate solutions ("WKB 
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solutions") that are usually numerically much less expensive than ab initio quantum 
mechanical calculations. In context of the propagator methods discussed here the semi- 
classical method can be described by the following steps: 

1. Find number of paths Adjoining source r' and destination r. This is an effective 
method for finding the caustics (turning surfaces). 

2. Find all trajectories r^t) leading from r' to r. 

3. Establish classical weight Pfc(r), i. e., the local density of trajectories, for each 
path. 

4. Determine its semiclassical phase "t^r) from the reduced action along the path. 

5. Create semiclassical approximation to the Green function G sc (r,r';E). 

6. Use uniform approximations to correct the divergence of the WKB solution near 
the caustics. 

Therefore, the semiclassical method requires knowledge of the classical trajectories 
rj.(f), their local density pt(r), and the corresponding action fields ^(r^jE). A 
semiclassical treatment of the photodetachment problem (Sec. 13.4b can be found e. g. in 
Ref. Il26l . Recently, the technique has been applied to the more complicated dynamics 
of electrons in parallel electric and magnetic fields ll50ll5"Tl . Here, we concentrate on 
the related problem of semiclassical motion in the Hall configuration (Sec. |6]>. 

Classical orbits for the Hall effect: The propagation of electrons in crossed elec- 
tric and magnetic fields has received much attention in classical physics. Indeed, the 
corresponding trajectory field of electrons emitted by a point source located at the ori- 
gin (for convenience), plotted in Fig. Q3] looks interesting by itself. Again assuming 
the Hall geometry displayed in Fig. Q~2] the motion of the charges is governed by the 
Hamiltonian d43l . and we find the family of orbits: 

r(t) = ( VDt \ + — ( ~ Voy v 0x -v D \ (cos(2co L t)\ 1 ( voy \ (54) 
w V / 2c °l Vox - v D voy J \sm(2co L t)J 2co L \v D -v 0x ) ' y 

Here, the initial velocity vector depends on the emission angle 9 and is given by 

r(r = O) = vo=f Vo ^f VoCOS fl V (55) 

It must fulfil E = \mv\. Eq. ( |54l is conveniently interpreted as a sum of three terms: 
On average, only the first term contributes to the transport of the electron. The corre- 
sponding drift velocity (averaged over one cyclotron period T = Tt/(0£) reads 

vd=1 f +1 d/'r(f') = {£ X 3!)/<% 2 . (56) 
1 Jt 
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Figure 15: Caustics and the trajectory field (cf. Fig.Qj): The trajectory field traces out 
the caustics. Electrons pass through focal points. The number of paths connecting the 
source with a given point increases with the magnetic field. 



In a "drift" reference frame that is moving with this velocity the otherwise trochoidal 
orbit becomes a circle with angle-dependent radius R(9) whose center is shifted from 
the origin by a constant displacement r ( 

fi (e) 2 = vg-2vov fl cos9+vf, J_/ vorinfl \ 

y ' 4col 2co L \ v D -v cosO) 

Some sample orbits are plotted in Fig. Qj] Variation of the angle yields the trajectory 
field displayed in Fig. [TBI 



The Hall conductivity tensor: Noting that the velocity in turn is related to the clas- 
sical current density j 

S=Nev D , (58) 

where N denotes the electron density, one can extract the resistivity tensor p (or its 
inverse, the conductivity tensor a) from Ohm's law in the two-dimensional (x,y)- 
plane: 

J-p-v - --«-'-£(? o 1 ) 

This remarkable equation predicts a finite conductivity even in the absence of scatter- 
ing, which is usually invoked in theories of conduction in order to guarantee a finite 
carrier velocity. We note that the drift velocity is independent of the kinetic energy of 
the electrons. 



Closed orbits in the classical picture: Eq. d58l i shows that the current in the Hall 
conductor is determined by the density of states N. According to Sec. [3] the local 
density of states (LDOS) n(E) within a narrow energy interval is related to the Green 
function G(o,o;£) (f4TT >. which, from a semiclassical perspective, is governed by those 
trajectories that return to the source (r = r' = o). This is the basic motivation of the 
closed orbit theory [52 53] for source processes. 

These trajectories, which lead from the origin back to the origin, are best found by 
using the classical action. Here, it is convenient to start with the time-dependent action 
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Figure 16: Density of states in crossed electric and magnetic fields. The curves show 
the quantum (solid) and semiclassical result (dashed line). The staircase structure de- 
notes the count of closed orbits. Crossed electric field S = 4000 V/m and magnetic 
field 38 = 5 T. 



functional 5 c i(r,f |r',0). In crossed electric and magnetic fields, this classical action is 
uniquely given by 

S A {o,t\o,Q) = --v 2 D t + —cot{o) L t)v 2 D t 2 . (60) 

This expression describes the single closed orbit returning to the source in a prede- 
termined time of flight t. However, we are rather interested in the energy E of the 
electron, 

M d (or|o,0) 
at 

For fixed emission energy E, this is an implicit equation for the time of flight t, and 
generally several solutions t k , pertaining to distinct classical trajectories r^(t), exist. 
(To find their initial velocities vo, it is sufficient to set r = o and t = t k in the equation 
of motion ( 154b . and solve the ensuing linear equation system for vq x and vo v .) The 
reduced action for each contributing path then follows from the Legendre transform: 

W k (o,o;E) = S d (o,t k \o,0)+Et k . (62) 

As shown in Fig. [161 the number of closed orbits increases with the magnetic field 
strength. 

Density of states and propagator: The classical action d60b is an important ingredi- 
ent of the quantum-mechanical time-evolution operator. Using Eqs. (HTt and d42l . it is 
possible to relate the local density of states with the propagator via [40, 52ll54l 

4%( E ) = ^Th jyTe> ET / r 'K(o,T\o,0), (63) 
where the time-dependent propagator ( l44b is given by 

^°' r l°' ) = 9^^, x expfa c iMo,0)|. (64) 
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Figure 17: Principal structure of the classical action in the complex time plane for 
N = 3 saddle points (closed orbits). The dashed line denotes the integration path. 
Singularities are denoted by x and saddle points by )(. Note that the singularity at the 
origin arises from the prefactor in the propagator and not from the classical action at 
f = 0. 



As we have shown before, this expression can be evaluated analytically in terms of 
harmonic oscillator eigenstates in energy space (07]). 

Quantum result versus semiclassical approach: An asymptotic evaluation of the 
integral d63l ) provides the link between closed orbits and the density of states. The 
original path of integration follows the real time-axis. Analytic continuation of the 
propagator makes it possible to deform this path of integration to the one sketched in 
Fig- 03 This path passes through saddle points of the exponent (denoted by )( in the 
figure) using the paths of steepest descent. The singularities in the integrand at times 
T = kn/ (Ol (denoted by x) are avoided. The only contribution of a singularity comes 
from t = 0, which may be evaluated by the residue theorem: 

/origin = r^-r / df t iEt ' h K{o,t\o,Q) = (65) 

2k h J 2n % 

Comparison of Eqs. d63t and d64l l with ( t6Tb shows that the saddle points of the inte- 
grand coincide with the classical times of flight for the various closed orbits. Adding 
their contributions yields the semiclassical result: 

m(0 L N e iW k (o.o;E)/ri+ijis g n[S c ,(o,t k \o.0)}/4 

n T / / ■■ = (66) 

47T 2 i h- t\ sin(o) L r A .) A /|S c /(o,^|o,0)|/(2wft)_ 

Fig. Q~6] compares the semiclassical and quantum results. Despite their very different 
origins (sum over Landau levels vs. interfering classical trajectories), they are in strik- 
ing agreement. In the "plateau regions" of the conductivity, destructive interference 
between the properly weighted classical trajectories strongly suppresses the LDOS, 
leading to a quantization into separated levels with a substructure. Note that the num- 
ber of trajectories does not change in each Landau level. It is the relative phase that 
modulates the LDOS. 



„( 2D ) -/ . . +2 Re 

n sc,Sy*3$ ~ 'origin t^m 
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